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ABSTRACT 

— j . We study by Monte Carlo computer simulations random sequential adsorption 

c : 

O . (RSA) with diffusional relaxation, of lattice hard squares in two dimensions. While 

CJ ■ 

for RSA without diffusion the coverage approaches its maximum jamming value (large- 
rS I time fractional coverage) exponentially, added diffusion allows the deposition process 

to proceed to the full coverage. The approach to the full coverage is consistent with 
the ~ t~^/^ power law reminiscent of the equilibrium cluster coarsening in models with 
nonconserved order-parameter dynamics. 



PACS numbers: 68.10.Jy, 02.50.+S, 82.65.-i 



Random sequential adsorption (RSA) models have been studied extensively due to 
their relevance to deposition processes on surfaces [1]. The depositing particles are rep- 
resented by hard-core extended objects; they are not allowed to overlap. In monolayer 
deposition of colloidal particles and macromolecules [2] one can further assume that the 
adhesion process is irreversible. Once a particle is in place its relaxation on the surface 
proceeds on time scales much larger than the deposition process. 

However, recent experiments on protein adhesion at surfaces [3] indicate that in 
these systems effects of surface relaxation, presumably due to diffusional rearrange- 
ment of particles, are observable on time scales of the deposition process. The resulting 
large-time coverage is denser than in fully irreversible RSA and in fact it is experimen- 
tally comparable to the fully packed (i.e., locally semi-crystalline) particle arrangement. 
Studies of RSA with diffusional relaxation by analytical means encounter several diffi- 
culties associated with possible collective effects in hard-core particle systems at high 
densities (such as, for instance, phase separation), and with the possibility, in certain 
lattice models, of locally "gridlocked" vacant sites. The latter effect may actually pre- 
vent full coverage in some models; this matter remains an open problem at this time. 

Both difficulties are not present in ID: there are no equilibrium phase transitions, 
traces of which might manifest themselves as collective effects in D > 1 deposition 
with diffusion, and furthermore diffusional relaxation leads to simple hopping-diffusion 
interpretation of the motion of vacant sites in ID which recombine to form larger open 
voids accessible to deposition attempts. Thus, both extensive numerical studies and 
their analytical interpretation were possible in ID [4]. For higher-D models however, 
no results were reported in the literature, to our knowledge. 

In this work we report extensive numerical simulations of the RSA process with dif- 
fusional relaxation, for the lattice hard-square model [5], i.e., the square-lattice hard-core 
model with nearest-neighbor exclusion. This model is well studied for its equilibrium 



phase transition [5] which is second-order with disordered phase at low densities and 
two coexisting ordered phases, corresponding to two different sublattice particle cov- 
erage arrangements, at high densities. Another simplifying feature of the hard-square 
model is that the only possible gridlocked (locally frozen) vacancies are parts of domain 
walls (see further below). As a result the coverage reaches the full crystalline limit at 
large times, by a process of diffusional domain wall motion leading to cluster growth 
reminiscent of quenched binary alloys and fluids at low temperatures [6] . 

The approach to the full coverage is numerically consistent with the ~ t~^/^ law 
which in turn can be related to the domain size growth with time as t^'^ as expected 
for order-parameter nonconserving dynamics. In the remainder of this work we first 
report computational details and numerical results. We then describe the domain- 
growth interpretation of the dynamics at high densities, as well as discuss some other 
collective effects observed in our simulations. 

In each Monte Carlo trial of our simulation on a L x L square lattice with periodic 
boundary conditions, a site is chosen at random. Then with probability r we attempt to 
deposit a particle and with probability (1 — r) we try diffusion. In the case of deposition, 
we check if the chosen site and its four nearest-neighbor sites are all empty. If indeed 
they are empty the deposition is performed. The chosen site is marked as occupied. 
If however any of the five sites are already occupied, then the deposition attempt is 
rejected and the configuration remains unchanged. In the case of diffusion, which is of 
course possible only if there is a particle at the selected site, we choose at random with 
equal probability a direction (up, down, left, or right) and try to move this particle 
by one lattice spacing. A move is made if the targeted new site and its three nearest 
neighbors different from the "source" site, are all empty. If the attempted move is not 
possible, the particle stays at its original position. 

A unit Monte Carlo time step is defined such that each lattice site is checked once 



on average. This corresponds to L? trials as described earlier. This time scale, T, is 
conveniently related to the physically more interesting time t defined to have fixed de- 
position attempt rate per site, with varying relative diffusion attempt rate proportional 

to (1 — r) /r, 

t = rT. (1) 

We implement the dynamics in two different ways. The first one is a straightforward 
simulation of the model. However, it is inefficient at late stages, when most trials are 
rejected. Thus we have implemented the same dynamics with an event-driven method 
[7]. We keep a list of all the possible moves, deposition and diffusion. Then we pick a 
move according to proper probability and always carry it out. The list is then updated if 
necessary. The Monte Carlo time is incremented according to AT = —{lnx)/R, where 
X is a uniform random number between and 1, while R is the rate of the system 
configuration changes in the original dynamics. 

Suppose our list contains A^depo sites available for deposition and A^diff possible 
diffusional moves. Each occupied lattice site contributes 0, 1, 2, 3 or 4 possible diffusional 
moves, while each empty lattice site contributes or 1 deposition counts, all depending 
on the nearest-neighbor site configuration. Then R is calculated as 

R = rNaepo + \{l-r)Nam. (2) 

Numerical estimates were obtained for the following quantities. The coverage, 0{T), 
was defined as the total number of sites occupied divided by L^/2. The deposition 
process always began with empty substrate so that the coverage 9 increased from to 
1 at full saturation. The "susceptibility" measuring fiuctuations of the magnetization, 
was defined by 

X = L'[{m')-{\m\f], (3) 
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where the average ( ) is over independent runs. The magnetization or order parameter 
was defined as usual [5] by assigning "spin" values +1 to particles on one of the sublat- 
tices and —1 on another sublattice. Empty sites on both sublattices were not counted 
(eflFectively having spin values 0). Thus the magnetization magnitude \m\ is defined as 
the difference of the number of particles deposited on two sublattices, normalized by 
the lattice size L^. The values of |rM| are thus from to ^. 

The effective domain size, i{T), was defined as in equilibrium-model studies of 
cluster coarsening [6], by 



e = 2Lv/(m2) . (4) 

Here the normalization is such that a uniform state (single-domain) gives the size of the 
system, L. 

Series of snapshots of the coverage buildup are shown in Fig. 1. As is usually done 
for the equilibrium hard-square system [5] , particles are represented by squares of size 
v2 X v2, rotated 45° with respect to the original square lattice on which the particle 
centers are deposited. Particles on the even and odd sublattices (the sum of the x and 
y coordinates even or odd) are shown in different shades. 

The time-dependence of the coverage is illustrated in Fig. 2. The general features 
are similar to those found in the ID studies [4] . For fixed deposition rate, corresponding 
to the time scale t defined in (1), added diffusion [rate ~ (1 — r)/r] always speeds up 
coverage growth. For high coverage, the r < 1 plots are reminiscent of domain growth 
in phase-separation models [6]. 

This similarity with domain growth dynamics can be made more quantitative. Let 
us first consider the coverage 0{T) for large times. In the case of deposition rate r = 1 
(no diffusion), the approach to the jamming coverage ^(oo) ~ 0.728 < 1 [8] is generally 
exponentially fast for lattice models [9]. With diffusion, one can always reach the full 
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coverage 9{oo) = 1. However, the approach to the fuU coverage is slow, power-law, as 
indicated in Fig. 2. In ID, the power-law behavior was related to the coverage growth at 
large times by the process of hopping and recombination (opening up deposition sites) 
of small empty regions [4] . 

The coverage growth mechanism for large times, in the 2D hard-square model is 
instead due to interfacial dynamics. As illustrated in Fig. 1, the void space at late 
times consists of domain walls separating spin-up and spin-down ordered regions. Since 
a typical domain has area ~ £^(T) and boundary ~ i{T), we anticipate that for large 
times 

1 - 9{T) (X r\T) . (5) 

Indeed, 1 — 9 is just the void area fraction. The large-time behavior of the domain size 
in models with nonconserved order-parameter dynamics is typically diffusional ~ T^/^. 
For £{T) we report the direct numerical verification later. For the coverage, we found 
that the data roughly fit the power law, 

1 - 9{T) oc T-i/2 , (6) 

for T > 10^ in typical runs such as illustrated in Fig. 2. Thus, the RSA quantity 1 — 9(T) 
behaves analogously to the energy excess in equilibrium domain growth problems. 

However, on careful examination we noted that the slopes obtained by least-square 
fits of the straight-line portions of the curves for large T were not exactly —0.5 but rather 
ranged from about —0.48 to about —0.52 when the rate parameter r was decreased from 
0.8 to 0.01, corresponding to accelerating the relative diffusion rate. These differences 
could not be fully attributed to statistical errors. Presumably, they represent effects of 
corrections to the leading power law behavior. 

The "susceptibility" x for a given finite size L has a peak and then decreases to 
zero, indicating long-range order for large T; see Fig. 3. The peak location seems size- 
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dependent, at Tpeak oc L^; thus, it is difficult to observe for large system sizes. Since 
finite-size effects set in for i{T) ~ L, which given the "bulk" power law i{T) ~ T^/"^ 
leads precisely to the criterion T r^ L"^, we expect this maximum in fluctuations to be a 
manifestation of the ordering process at high densities. 

For equilibrium hard squares [5] , the critical-point density corresponds to coverage 
9 ~ 0.742. However, in all our simulations, even with the fastest relative diffusion rate 
(the case r = 0.01), we found no interesting features in x(T) for times for which the 
coverage was near the critical value. Since numerical effort to reach a given coverage 
increases with increasing the diffusion rate at the expense of deposition attempts, there 
still remains a numerical challenge to observe the buildup of rounded critical-point 
fluctuations in x or other quantity, for RSA with diffusion. Our present simulations turn 
out to be sensitive only to those collective effects which are associated with ordering at 
higher than critical-point coverages. 

In Fig. 4, the effective domain size i{T) is plotted vs. time for r = 0.1. The 
large-time asymptotic law is well fitted by T^/^. As already mentioned, the late-stage 
dynamics of our model is similar to the kinetics of ordering by quenching an equilibrium 
system into a two-phase region. In fact, our model can be mapped approximately to the 
Ising model at zero temperature. Note that the straight-line, diagonal sections of the 
domain boundaries separating two phases, can not move. Only the corners can evolve 
due to diffusion, just like corners of interfaces in the Ising models at zero temperature 
evolve by order-parameter nonconserving dynamics. It has been known that the t^''^ 
law is robust, independent of the number of ordered phases, dimensionality, and details 
of interactions in domain growth problems [6] . 

In summary, we found that when diffusion is introduced in random sequential ad- 
sorption processes, the system can relax to the full coverage, consistent with experiments 
on protein adsorption at surfaces. The approach to the full coverage is slow (power law). 



The dynamics for large coverages is governed by domain growth. 
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FIGURE CAPTIONS 



Fig. 1. Configurations of ttie tiard-square RSA-witii-diffusion model obtained with 
the deposition rate parameter r = 0.1, at times T = 10, 20, 30, and 100. Particles 
centered on the even sublattice are shown in grey, while those on the odd sublattice 
are shown in black. Lattice size was 32 x 32. 

Fig. 2. Void area fraction, 1 — 0, plotted against time t = rT on a double-logarithmic 
scale, for the deposition rate parameter r = 1, 0.8, 0.1, 0.01, and system size 
521x521. The data were averaged over 800, 600, 800, and 400 runs for r = 1, 0.8, 
0.1, and 0.01, respectively. The dashed-dotted line has slope of — |. 

Fig. 3. Fluctuation of the magnetization, %, vs. time T, for system sizes L =4, 8, 16, 
32, and 64. The deposition rate was r = 0.1. The data were averaged over 10^ to 
10^ runs. 

Fig. 4. Effective domain size i, plotted vs. time T on a double-logarithmic scale, with 
r = 0.1 and system size 256x256. The data were averaged over 1500 runs. 
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